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ABSTRACT 

The high angular resolution and sensitivity of the Chandra X-ray Observatory has yielded large num- 
bers of faint X-ray sources with measured redshifts in the soft (0.5 — 2 keV) and hard (2 — 8 keV) energy 
bands. Many of these sources show few obvious optical signatures of active galactic nuclei (AGN). We 
use Chandra observations of the Hubble Deep Field North region, A370, and the Hawaii Survey Fields 
SSA13 and SSA22, together with the ROSAT Ultra Deep Survey soft sample and the ASCA Large Sky 
Survey hard sample, to construct rest-frame 2 — 8 keV luminosity functions versus redshift for all the 
X-ray sources, regardless of their optical AGN characteristics. At z = 0.1 — 1 most of the 2 — 8 keV light 
density arises in sources with luminosities in the 10 42 ergs s' 1 to 10 44 ergs s _1 range. We show that 
the number density of sources in this luminosity range is rising, or is at least constant, with decreasing 
redshift. Broad-line AGN are the dominant population at higher luminosities, and these sources show 
the well-known rapid positive evolution with increasing redshift to z ~ 3. We argue that the dominant 
supermassive black hole formation has occurred at recent times in objects with low accretion mass flow 
rates rather than at earlier times in more X-ray luminous objects with high accretion mass flow rates. 

Subject headings: cosmology: observations — galaxies: evolution — galaxies: formation — galaxies: 
active 



1. INTRODUCTION 

Hard X-ray surveys provide the most direct probe of 
■ supermassive black hole (SMBH) accretion activity, since 
high energy X-rays can penetrate extremely large column 
densities of gas and dust. Early analyses of deep Chandra 
X-ray surveys found that optical selections of active galac- 
tic nuclei (AGN) had missed large numbers of accreting 
SMBHs. Barger et al. (2001b) showed that the redshift 
history of SMBH growth is more strongly peaked to low 
redshifts than had been inferred from optically-selected 
samples. 

With the high spatial resolution and energy sensitiv- 
ity of deep Chandra observations, X-ray samples can be 
selected at both low and high redshifts in the same rest- 
frame hard energy band and their optical counterparts, 
and hence spectroscopic redshifts (~ 65% can be identi- 
fied), unambiguously determined. The hardest band for 
which the hard X-ray luminosity function (HXLF) can 
presently be determined over a wide range of redshifts 
is rest-frame 2 — 8 keV. Although some Compton-thick 
sources may still be omitted, the accounting will be far 
more complete than that of any other available sample. 

In this paper we use the 1 Ms exposure of the Chandra 
Deep Field North (CDF-N) and three - 100 Ks expo- 
sures of the A370, SSA13, and SSA22 fields, together with 
previous ROSAT and ASCA data, to construct the rest- 
frame 2 — 8 keV luminosity function from z = to 4 and 
model the SMBH accretion history. We assume Hm = 1/3, 
ft A = 2/3, and H Q = 65 km s' 1 Mpc' 1 . We use L x to 
denote the rest-frame 2 — 8 keV luminosity. 



2. X-RAY SAMPLE SELECTION 

The faintest data are taken from the 1 Ms CDF- 
N catalogs of Brandt et al. (2001; hereafter B01) and 
Barger et al. (2002). We exclude known stars and the 
small number of extended X-ray sources. The flux lim- 
its are « 3 x 10' 17 ergs cm' 2 s' 1 (0.5 — 2 keV) and 
w 2 x 10' 16 ergs cm' 2 s' 1 (2-8 keV). Below a flux 
of 1.6 x 10' 14 ergs cm' 2 s' 1 (2 — 8 keV) and 6 x 
10 -15 ergs cm' 2 s' 1 (0.5 — 2 keV) we include only sources 
lying within 8' of the aimpoint where the redshift iden- 
tifications are most complete; above these fluxes we in- 
clude all the sources. This restricted sample contains 199 
(52%) sources in the 2 - 8 keV band and 219 (57%) in 
the 0.5 — 2 keV band, where the parentheses are the per- 
centages with redshifts. We supplement the sample at 
fluxes above 2 x 10 -15 ergs cm' 2 s' 1 (2 — 8 keV) and 
3 x 10' 16 ergs cm" 2 s" 1 (0.5 — 2 keV) with data from the 
SSA13 field (Barger et al. 2001a), the A370 field (Barger 
et al. 2001b), and the SSA22 field (A. J. Barger et al., in 
preparation). This adds a further 67 (72%) sources in the 
2 - 8 keV band and 88 (60%) in the 0.5 - 2 keV band. We 
further supplement the sample with the larger area bright 
flux samples in the soft band from ROSAT (Lehmann et 
al. 2001) and the hard band from ASCA (Akiyama et al. 
2000). These have highly complete spectroscopic identifi- 
cations. 

For the conversion of counts to flux in the CDF-N, we 
use the X-ray spectral photon indices T determined in B01 
from the hard-band to soft-band counts ratios. For the 
conversions in the remaining Chandra fields, we assume a 
fixed r = 1.2, which is approximately the average T in 
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the B01 sample. For the flux conversions in the ROSAT 
and ASCA samples, we use the softer T values assumed by 
Lehmann et al. (2001) and Akiyama et al. (2000), which 
are appropriate for higher flux sources. 

We determined the solid angle covered by the sample 
at a given flux by comparing the observed numbers of 
objects as a function of flux with the measured number 
counts in the appropriate energy band. For 2 — 8 keV we 
use the averaged counts given in Cowie et al. (2002), and 
for 0.5 — 2 keV we use the power-law fits of Mushotzky et 
al. (2000), which agree with other determinations (e.g., 
BOf; Rosati et al. 2002). The solid angle covered by 
the combined 2 — 8 keV samples ranges from just un- 
der 0.01 deg 2 at the faintest fluxes to 5.8 deg 2 at the 
highest fluxes. At 10~ 14 ergs cm~ 2 s _1 (2 — 8 keV) the 
solid angle is 0.11 deg 2 . The maximum solid angle cov- 
ered by the combined 0.5 — 2 keV samples is 0.42 deg 2 . 
At 10~ 15 ergs cm~ 2 s _1 (0.5 - 2 keV) the solid angle is 
0.12 deg 2 . 




Fig. 1. — Redshift versus flux for sources selected in (a) the 
0.5 — 2 keV band and (b) the 2 — 8 keV band. Solid curves show flux 
versus redshift for sources with L x = 10 42 ergs s _1 (lowest curve), 
10 43 ergs s — 1 , 10 44 ergs s _1 , and 10 45 ergs s~ 1 (highest curve), 
computed with a _R"-correction for a T = 1.8 power-law spectrum. 
Dashed lines show the redshift intervals used in constructing the 
HXLFs. Sources with broad-line optical spectra are enclosed in a 
second, larger symbol. Spectroscopically unidentified sources are 
denoted by plus signs at z = —0.3. 

In Fig. la and b we show redshift versus flux for 
the soft and hard selected samples, respectively, includ- 
ing sources that have not been spectroscopically identi- 
fied. The solid curves correspond to loci of constant L x . 
Sources with broad- line optical/UV spectra are denoted 



by a second, larger symbol. Any source more luminous 
than 10 42 ergs s" 1 is very likely an AGN on energetic 
grounds (Zezas, Georgantopoulos, & Ward 1998; Moran 
et al. 1999), though many of the intermediate luminosity 
sources do not show obvious AGN signatures in the optical. 
It is clear from Fig. 1 that the highest X-ray luminosity 
sources are mostly broad-line AGN. 

3. EVOLUTION OF THE AGN X-RAY LUMINOSITY 
FUNCTION 

In the present section we concentrate on two redshift 
intervals: z = 0.1 — 1, as shown by the dashed lines in 
Fig. lb, and z = 2 — 4, as shown in Fig. la. With these 
redshift intervals we can compare samples chosen in ap- 
proximately the same rest-frame hard X-ray energy range 
of 3 — 12 keV (the observed 2 — 8 kev band at the center 
of the low redshift interval) and 2 — 8 keV (the observed 
0.5 — 2 keV band at the center of the high redshift in- 
terval). The remaining small differential if-corrections to 
rest-frame 2 — 8 keV are made by assuming that the spec- 
tra can be approximated by T = 1.8 power-laws. We have 
chosen to use this spectral index since the differential K 
correction is mostly significant in correcting from higher 
to lower energies where the spectra may be better repre- 
sented by the unabsorbed intrinsic power law (see Barger 
et al. 2002). However, the results are only weakly sensitive 
to this assumption. For aT = 1.2 power-law K-correction 
the HXLF at z = 2 — 4 is essentially unchanged while the 
HXLFs at z = 0.1 — 1 shifts to slightly lower luminosities 
(a factor of 1.27) while the shape remains approximately 
the same. 

We have chosen not to correct for the intrinsic absorp- 
tion since this too requires an assumption about the in- 
trinsic spectrum and also about the form of the absorp- 
tion and any scattering. Since we are comparing in the 
same rest frame at different redshifts this will not affect 
the HXLF comparison unless there is substantial evolution 
in the intrinsic absorption. Moreover in the 2 — 8 keV rest 
frame band the average absorption corrections are not very 
large. If we assume an intrinsic spectral index of T = 1.8 
and a simple photoelectric absorption the average absorp- 
tion correction to L x is 1.14 in the 0.1 — 1 redshift range 
and 1.53 in the z = 2 — 4 range. 

We probe to an L x below 10 43 ergs s" 1 in the high red- 
shift interval and below 10 42 ergs s _1 in the low redshift in- 
terval. Both samples are highly complete in redshift iden- 
tifications at high luminosities. At L x > 10 44 ergs s _1 in 
the high redshift interval there are 22 spectroscopically 
identified sources, 21 of which are broad-line AGN. At 
most there are a further 18 spectroscopically unidentified 
sources that could lie at these luminosities in this red- 
shift interval; however, a substantial fraction of these are 
likely at other redshifts. At L x > 10 44 ergs s _1 in the low 
redshift interval, 85% of the identified sources are broad- 
line AGN, and there are only three unidentified sources 
that could lie at these luminosities in this redshift inter- 
val. Thus, most of the high luminosity X-ray sources are 
broad-line AGN rather than obscured AGN. 

The differential HXLF as a function of X-ray luminosity 
and redshift, d$(L x , z)/d\ogL x , is the number of X-ray 
sources per unit comoving volume per unit logarithmic 
luminosity. We adopt the 1/V a method to estimate the 
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binned HXLF in different redshift shelis as 
d log Ac AlogL x ' 



(1) 



where the is the total comoving volume in the redshift 
bin Az in which source i with luminosity L x could lie and 
still be included in the sample and the summation is over 
all sources in the given redshift interval and luminosity 
bin. (Log denotes the base 10 logarithm.) 
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Fig. 2. — Rest-frame 2 — 8 keV luminosity function per unit loga- 
rithmic luminosity in the redshift intervals z = 0.1 — 1 (open squares) 
and z = 2 — 4 (solid diamonds). In the z = 0.1 — 1 (z = 2 — 4) 
interval the HXLF was computed from observed-frame 2 — 8 keV 
(0.5 — 2 keV). An intrinsic T = 1.8 was assumed, for which there 
is only a small differential A"-correction to rcst-framc 2 — 8 keV. 
Poissonian lcr uncertainties are based on the number of sources in 
each luminosity bin. Dotted (z = 0.1 — 1) and dashed (z = 2 — 4) 
lines arc maximal HXLFs found by assigning all the spectroscop- 
ically unidentified sources to the center of each redshift interval. 
Dot-dashed line is an upper bound on the z = 2 — 4 HXLF ob- 
tained by selecting in observed-frame 2 — 8 keV and then assigning 
all the unidentified sources to this redshift interval. Solid lines show 
a broken power-law fit to the low redshift HXLF with a slope of 
—0.6 below 1.2 X 10 43 ergs s — 1 and —1.35 above, and to the high 
redshift HXLF with a slope of 0.0 below 4.8 X 10 44 ergs s _1 and 
—2.35 above. The upper L x fit in the high redshift interval is based 
on the HXLF of Miyaji ct al. (2001) shown in Fig. 3. 

Our measured HXLFs for the two redshift intervals are 
shown in Fig. 2 with Poissonian lcr uncertainties based on 
the number of galaxies in each luminosity bin. However, 
since incompleteness is a potentially larger source of error, 
we recomputed the HXLFs by assigning all the spectro- 
scopically unidentified sources to the center of each red- 
shift interval (i.e., we included all the unidentified sources 
in both redshift intervals). Dotted (z = 0.1 — 1) and dashed 
(z = 2 — 4) lines in Fig. 2 show these maximal limits. Be- 
cause the spectroscopic identifications are much more com- 
plete at the high X-ray fluxes, the associated systematic 
uncertainties are larger at lower L x . In principle photo- 
metric redshifts can be estimated for the sources, and at 
higher redshifts this procedure is effective (A. J. Barger 
et al., in preparation). However, over the z = — 4 range 
there could be considerable uncertainties in assigning pho- 
tometric redshifts to sources with substantial AGN con- 
tributions, so we prefer our more conservative approach. 
Analyses of the HXLF using alternate techniques suggest 
that the lowest luminosity point in the determinations may 
be biased low (e.g. Page and Carrera 2000, Miyaji et al. 
2001) but the effect is small compared to the statistical 
and systematic uncertainties. 



At L x > 10 ergs s~ 4 Fig. 2 shows the familiar rise 
of the number counts with redshift. These counts are 
dominated by broad-line AGN. However, at lower L x the 
z = 0.1 — 1 counts lie above the z = 2 — 4 counts, or 
are comparable if we use the maximal correction for in- 
completeness. Since it is extremely unlikely that all the 
unidentified sources lie in the z = 2 — 4 range (this would 
result in a large deficiency at z = 1 — 2), the present-day 
number density of sources at L x ~ 10 43 ergs s -1 is higher 
than it was in the past. This result cannot be significantly 
changed by highly obscured sources: the dot-dashed line 
shows the upper bound on the z — 2 — 4 HXLF if we 
selected in observed-frame 2 — 8 keV, which corresponds 
to rest-frame 8 — 32 keV, and placed all the unidentified 
sources in this redshift interval. The resulting upper limit 
is a factor of ~ 2 higher than the upper limit obtained 
from the observed-frame 0.5 — 2 keV selection. 

All the X-ray sources are included in Fig. 2 without 
regard to optical spectroscopic classification. Previous 
HXLF analyses have generally focused on broad-line AGN 
where the counterparts can be determined even in low res- 
olution X-ray data (e.g., Boyle et al. 1998). To make a di- 
rect comparison with these earlier studies, in Fig. 3 we plot 
our spectroscopically identified broad-line AGN HXLFs in 
the intervals z = 0.1 — 1 and z = 2—4 with those derived by 
La Franca et al. (2002; z = 0.5) from BeppoSAX data and 
by Miyaji et al. (2001; z = 2 -4) from ROSAT data. Our 
HXLFs agree with these brighter HXLFs in the overlap re- 
gion where there is also some overlap in the datasets used. 
At L x ~ 10 43 ergs s _1 there is little apparent evolution in 
the broad-line AGN HXLF, while at higher luminosities 
the HXLF increases with redshift. 
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Fig. 3. — Rest-frame 2—8 keV broad-line AGN luminosity function 
per unit logarithmic luminosity in the redshift intervals z = 0.1 — 1 
(solid squares) and z = 2 — 4 (solid diamonds). Dashed line is the 
maximal HXLF in the z = 2 — 4 interval found by assigning all 
the spectroscopically unidentified sources to the center of the in- 
terval. However, we expect that most of the broad-line AGN have 
been identified since they are optically bright and straightforward 
to identify spectroscopically. Open boxes (open diamonds) show the 
broad-line AGN HXLF at z = 0.5 (z = 2 - 4) derived by La Franca 
et al. (2002) (Miyaji et al. 2001). 



4. COSMIC EVOLUTION 

In Fig. 4a we show the comoving number density evo- 
lution with redshift of L x > 10 42 ergs s _1 sources. At 
these luminosities the probability that the X-ray sources 
are powered by anything other than AGN is extremely 
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low, so we are mapping the evolution of intermediate lu- 
minosity AGN. In the z = 1 — 2 rcdshift bin we show 
the number densities from both observed-frame 2 — 8 keV 
(open squares) and observed-frame 0.5 — 2 keV (solid dia- 
monds); the results are similar, although the soft band may 
be missing some obscured sources. The solid (2 — 8 keV) 
and dotted (0.5 — 2 keV) horizontal bars show the number 
densities obtained if all the spectroscopically unidentified 
sources are assigned redshifts at the center of each redshift 
bin. The bars are not consistent with one another because 
all the unidentified sources are included in all the bins. 



We parameterize the evolution of X x as 
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Fig. 4. — Evolution with redshift of the comoving number density 
(a) and the 2 — 8 keV comoving luminosity density production rate 
of L x > 10 42 ergs s _1 sources (b). Open squares (solid diamonds) 
are measured from observed- frame 2 — 8 keV (0.5 — 2 keV). Pois- 
sonian la uncertainties are based on the number of sources in each 
bin. Dashed (0.5 — 2 keV) and solid (2 — 8 keV) bars show upper 
limits found by assigning all the unidentified sources to the center 
of each redshift bin, so the bars are not consistent with one another. 
To minimize this effect at high redshifts, we use a single 2 = 2 — 4 
redshift bin. At the highest redshifts a small correction (less than 
20%) has been applied to extrapolate from the minimum measured 
luminosity to L x = 10 42 ergs s — 1 . Note that the uncertainty in the 
\ x , which is weighted to higher luminosity sources, is smaller than 
the uncertainty in the number density. 

The number density of intermediate luminosity AGN 
does not rise rapidly from z ~ to z ~ 2.5 — 3, in con- 
trast to that of higher luminosity AGN (e.g., Miyaji et al. 
2000). In fact, the evolution with increasing rcdshift is flat 
or declining. Thus, AGN traced by X-ray luminosity show 
their own version of downsizing, just as star-forming galax- 
ies do (Cowie et al. 1996). High X-ray luminosity sources, 
which strongly overlap with broad- line AGN (quasar) pop- 
ulations, peak at higher redshifts, while intermediate lu- 
minosity sources are still common now. 

In Fig 4b we show the evolution with rcdshift of the 
2 — 8 keV comoving luminosity density production rate of 
L x > 10 42 ergs s _1 sources (X x ). Integration of the power- 
law fits to the HXLFs shown in Fig. 2 give a similar an- 
swer: 1.9 x 10 39 ergs s -1 Mpc~ 3 in the z = 0.1 — 1 redshift 
interval and 2.2 x 10 39 ergs s _1 Mpc~ 3 in the z = 1 — 4 in- 
terval. The largest uncertainty is the redshift distribution 
of the unidentified sources rather than the small effects of 
extrapolation outside our observed luminosity range. This 
incompleteness uncertainty is less than a factor of 2.9 even 
in the high redshift interval (Fig. 4b). 



X x 



1 + z 
1.5 



(2) 



where a ranges from —0.4 (all unidentified sources in 
the low redshift interval) to 1.1 (all unidentified sources 
in the high redshift interval). If only identified sources 
are included, \ x is roughly flat until the lowest z = 
0.1 — 0.5 interval in Fig. 4 where it may drop slightly and 
A « 2 x 10 39 ergs s _1 Mpc~ 3 . Thus, the 2-8 keV luminos- 
ity density production rate evolves slowly with redshift, as 
was first noted in Barger et al. (2001b). 

The SMBH accretion rate density is now given by 



Psmbh — B 



Ax(l-e) 



(3) 



where B is the correction from L x to the bolometric lu- 
minosity of the AGN, and e is the radiative efficiency of 
the accretion flow. Setting B ~ 40, which seems to 
be roughly valid for both unobscured and obscured AGN 
(Elvis et al. 1994; Barger et al. 2001b), and e ~ 0.1, which 
is close to the maximum possible radiative efficiency, 



Psmbh — C 



1 + z 
1.5 



(4) 



where C « 1.3 x 10~ 5 Mq yr _1 Mpc~ 3 . Even within 
the uncertainty in a, SMBH growth is continuing strongly 
to the present time. Integrating Eq. 4 for a — gives a 
present universal SMBH density of w 2 x 10 5 M Mpc -3 , 
comparable to the w 2 — 4 x 10 5 Mq Mpc -3 measured 
from local galaxies (Ferrarese & Merritt 2000; Gebhardt 
et al. 2000; Yu & Trcmaine 2002). The late formation of 
the X-ray background reduces the required local SBMH 
mass density of previous estimates that assumed the bulk 
of the production was earlier (Fabian and Iwasawa 1999). 

In summary, our results show that while the higher ac- 
cretion mass flow rates that power the most luminous AGN 
peaked at higher redshift and are now much rarer, lower 
accretion mass flow rates are still common at the present 
time. These less luminous events dominate SMBH for- 
mation. The longer duration of the lower mass flow rate 
systems provide constraints on theoretical modelling of the 
SMBH formation (e.g., Haehnelt, Natarajan, & Rees 1998; 
Haiman & Menou 2000; Yu & Tremaine 2002) and on the 
evolution of the accretion disks feeding the sources, which 
may evolve more rapidly if they are initially of higher mass 
(Duschl 2002, private communication). 
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